Quantum vortex tunneling in Y Ba2CuzO'j-s thin films 
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Cuprate films offer a unique opportunity to observe vortex tunneling effects, due to their unusu- 
ally low superfluid density and short coherence length. Here, we measure the magnetoresistance 
(MR) due to vortex motion of a long meander line of a superconducting film made of underdoped 
Y Ba2Cuz07-s. At low temperatures (T), the MR shows a significant deviation from Arrhenius 
activation. The data is consistent with two dimensional Variable Range Hopping (VRH) of single 
vortices, i.e. MR oc ea:p[— (Tq/T)^^^]. The VRH temperature scale To depends on the vortex tun- 
neling rates between pinning sites. We discuss its magnitude with respect to estimated parameters 
of the meander thin film. 
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Very soon after the discovery of the high temper- 
ature superconductors, the resistive transition was 
observed to broaden under magnetic fields, instead of 
shifting to lower temperatures This was explained 
as due to thermally activated vortex motion at high 
temperatures just below Tc, which gave rise to an 
induced voltage across the superconductor and thereby 
to a flux flow resistance Rff. At any given magnetic 
fleld in this regime, this resistance has the form of an 
Arrhenius law Rff (x exp[—{UQ/kBT)], where Uq is 
the activation energy which was discussed by several 
authors [1, H, Generally, flux flow and flux creep 
are possible at high temperatures where the pinning is 
relatively weak compared to the thermal energy. At low 
temperatures, where the pinning is stronger and thermal 
activation is much weaker, the dominant mechanism for 
flux motion is via quantum tunneling. Measurements 
of magnetic relaxation and transport by Stein et al. 
have shown a signature of quantum flux creep in 
Yi^xP'i'xBa2Cu30'j^s crystals. This behavior, where 
the creep is temperature independent, was observed in a 
very limited temperature range of about 2-4 K. Vortex 
tunneling in a 2D superconductor at temperatures much 
lower than the transition temperature Tc was discussed 
theoretically by Fisher, Tokuyasu and Young (FTY) 
0, and more recently by Auerbach, Arovas and Gosh 
(AAG) FTY studied quantum vortex tunneling in 
2D fllms at T close to zero near the superconductor 
to insulator glass transition. They found that this 
tunneling occurs via variable range hoping (VRH), but 
instead of the usual 1/3 power law in the exponent, at 
high flelds they predicted a different behavior where the 
exponent ranges between 2/3 and 4/5. This resulted 
from taking into account vortex-vortex interactions. At 
low flelds {H << Hc2), AAG calculated the tunneling 
rate of a single vortex between two pinning sites and 
the resulting flux tunneling resistivity. They found 
that this resistivity depends on temperature as the well 
known VRH in 2D, namely p cx exp[-{To/Ty/^]. In the 
present study we set up an experiment to test this AAG 
prediction in a thin fllm of Y Ba2CuzOT-s (YBCO) 
patterned into a very long and narrow meander line. We 
note that in a short microbridge of YBCO of typically 
a few hundred /im length, and under a magnetic field 



of several Tesla, the voltage induced by moving vortices 
becomes immeasurably small, and in practical terms, a 
critical current develops already at about 10-20 K below 
Tc- This prohibits measurements of Rff at low bias 
and low temperatures. In contrast, in a much longer 
meander line, the total DC component of the voltage 
generated by moving vortices is large, and a de-facto 
resistive state persists down to very low temperatures, 
thus enabling the vortex tunneling study. 



A question still arises as for why quantum vortex 
tunneling has not been unequivocally identified in low 
Tc superconductors? We believe that this is due to 
the fact that conventional BCS superconductors are 
generally three dimensional and characterized by a 
large coherence length and high superfluid density, all 
of which greatly inhibit vortex tunneling. In contrast, 
underdoped cuprate films, such as investigated here, 
offer a unique system of low superfluid density and 
quasi 2D superconductivity with very short coherence 
lengths. According to a recent theoretical work in 
such two dimensional "bosonic" superconductors, vortex 
tunneling can be manifested in terms of deviation from 
Arrhenius thermal activation of the magnetoresistance. 
Therefore, observing the (Tq/T)^/'^ exponent would 
be an exciting signature of the bosonic character of 
cuprate superconductivity as well as a conflrmation of 
short length-scale vortex tunneling effects. Concerning 
previous experimental results, we note that the phe- 
nomenon of vortex quantum creep has been investigated 
by both magnetization relaxation and transport in 
three dimensional cuprate crystals Unfortunately, 
although a temperature independent relaxation time and 
voltage drops have been observed at low temperatures, 
these experiments did not necessarily probe the low 
current regime where tunneling (rather than above 
barrier classical diffusion) is expected to dominate. 
The transport results explicitly show non-saturation at 
low currents which indicates that no vortex tunneling 
occurred in the low current, Ohmic regime. In our 
experiments, we study thin films with low current bias. 
This ensures that deviation from Arrhenius activation 
which we find is indeed related to vortex tunneling and 
not to classical above-barrier motion. 
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In the present study, we used epitaxial c-axis oriented 
films of YBCO, with the magnetic field applied parallel 
to this axis and perpendicular to the substrate wafer. 
The films were prepared by 355 nm laser ablation 
deposition on (100) SrTiO^ wafers of lOx 10 mm^ area. 
We investigated films of 50 and 100 nm thickness. In 
order to reduce their contact resistance, some films were 
coated with a 20 nm thick gold layer deposited at a 
high temperature (780 °C) in oxygen ambient. This 
gold overlayer consists of loosely connected ball-like 
grains as seen by atomic force microscope images. Its 
contribution to the resistance versus temperature of 
the bilayer was non-metallic and very small. This 
was verified in a control experiment in which another 
bilayer of 20 nm gold on 100 nm YBCO was prepared 
with the YBCO doped to have a very low Tc. Thus, 
its resistance could be measured down to very low 
temperatures. This bilayer was measured as deposited, 
without patterning. After measuring the resistance, the 
Au layer was removed using Ar ion milling, and the 
resistance was measured again. A simple calculation 
of two parallel resistors shows that the resistivity of 
the gold layer below 10 K was above 5 mil cm. This 
value is much higher than the resistivity of the 60 K 
YBCO phase which is about 0.1 mflcm at 100 K. 
Considering also the 1:5 thickness ratio between the 
gold and YBCO layers used, we find that the presence 
of the gold layer changed the resistance of the meander 
line by about 0.25%, which is negligible. Nevertheless, 
the aim of obtaining a low contact resistance with the 
YBCO film where the contact is made normal to the 
surface was certainly achieved. We note that the MR 
results measured with and without the gold layer were 
basically the same. However, those obtained with the 
gold overcoating were less noisy, and therefore we chose 
to present them here. 

The YBCO films and Au/YBCO bilayers were an- 
nealed in-situ in a controlled oxygen ambient to have 
a transition temperature of 50-60 K, in order to avoid 
apparent critical currents at low temperatures. We note 
that the presence of a critical current interferes with 
the resistance studies at low bias where the current 
should be proportional to the voltage. The films were 
patterned by deep UV lithography using a PMMA resist 
and Ar ion beam milling, into a 4 m long meander line. 
This meander line had 14 nm lines-width with 4 /xm 
lines-spacing on 8 x 9 mm? area of the wafer, while the 
remaining area was used for the four 1.5 x 1.5 mm^ 
contacts. In the patterning process, we limited the 
development time of the photo-resist in order to keep 
it continuous, but this left several shorts, effectively 
shortening the meander line. In addition, we found that 
in the patterning process a few defects were formed in 
the long meander line which we bridged with small silver 
paste dots of ~ 0.5 mm diameter. The gold overcoating 
layer was also very helpful in improving the contacts to 
these bridging dots. This bridging procedure and the 
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FIG. 1: (Color online) Resistance versus temperature of the 
YBCO meander line (of ~lm length and 100 nm x 14 /xm 
cross section area) under various cooling conditions with and 
without a magnetic field. The resistance of the meander line 
over a wider temperature range under 2 T field cooling is 
shown in the inset. 
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FIG. 2: (Color online) The resistance data of Fig. 1 plotted on 
a log scale versus inverse temperature. The standard vortex 
activation regions where R oc exp{ — Uo/kT) are shown for 
fields of 2 T and 6 T. 

shorts mentioned before led to a reduced effective length 
of our meander lines of about 1 m. This length was 
estimated from the measured resistance, the resistivity 
and the cross section area of the bare YBCO meander 
lines without the gold coating. For the transport mea- 
surements we used the standard four contacts technique 
in a dc mode, and cooling was done in a liquid helium 
cryostat with a base temperature of 2 K. 

The resistance versus temperature results of the 
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FIG. 3: (Color online) The magnetoresistance MR — R{H) — 
i?(0) derived from the data of Fig. 1, plotted on a log scale as 
a function of T~^^^ . The straight lines are linear fits to the 
data in the relevant low temperature regimes. In the inset the 

1/3 

slopes of these straight lines (on a In scale) Tq are plotted 
versus field, together with a fit to A.l/B^^^ . 



meander line are shown in Fig. 1 for zero field cooling 
(ZFC) and under magnetic fields of 1, 2, 4, and 6 T 
applied normal to the wafer. One can see the typical 
behavior characteristic of underdoped YBCO in the 
2-300 K temperature range in the inset to this figure. 
In the main panel, the superconducting transition 
temperature range is shown in more detail. The well 
known thermally activated broadening of the transition 
under increasing magnetic field is clearly observed, and 
in addition a significant resistive tail is seen already 
at zero field. To see the thermal activation functional 
dependence, and the data at low temperatures, the data 
of Fig. 1 was reploted in Fig. 2 on a log scale versus 
inverse temperature. One can see that the regimes 
where thermal activation holds are small, about 5 K at 
2 T and 7 K at 6 T. The onset of the superconducting 
transition under ZFC is T°"'"=* = 58iC, while that of 
the 6 T field cooled curve is shifted down to about 40 
K. This is qualitatively similar to the results of Palstra 
et ai, but with about twice the measured shift of the 90 
K phase of YBCO single crystals [1]. Also seen in Fig. 
2 is a knee in the ZFC curve at 42-34 K just below the 
main transition. This is a clear signature of weak links, 
which is very similar to the results on grain boundary 
junctions Q. These weak links can be due to defects in 
the STO substrate which are copied into the epitaxial 
film, presence of grains of the minority (~3%) a-axis 
oriented phase, lithography induced defects and so on. 
It is therefore clear that any possible vortex tunneling 
unmasked by other effects could be observed only at 
temperatures below about 20-30 K. 




FIG. 4: (Color online) The data of Fig. 3 plotted versus in- 
verse temperature to compare possible standard vortex acti- 
vation to vortex tunneling by VRH as in Fig. 3. The straight 
lines are guides to the eye. 

Next we study the temperature dependence of the mag- 
netoresistance of the meander line. We use the usual 
definition of MR = R{H) - R{0) where R{H) and R{0) 
are the resistances under a magnetic field H and at zero 
field, respectively. Using the magnetoresistance rather 
than the resistance under field itself has the advantage of 
eliminating all extrinsic effects such as those contributed 
by weak links and series resistors of the repairing silver 
paste dots. We can safely assume that these contribu- 
tions to the magnetoresistance are small since the rele- 
vant defects occupy a very small fraction of the area of 
the meander line, and their resistance has a much weaker 
dependence on the perpendicular magnetic field. In Fig. 
3 the MR is plotted versus T~^^^ where a linear depen- 
dence should indicate VRH in the YBCO planes, while in 
Fig. 4 the MR is plotted versus 1/T to test for a possible 
Arrhenius activation process. As one can see a linear be- 
havior is found in both figures at low temperatures, but 
the VRH range of Fig. 3 extends over 2-10 K which is 
about three times as large as the activation range of 2-5 
K as seen in Fig. 4. Although this temperature range is 
relatively small, the fact that the log{MR) cx T~^/^ be- 
havior exists in a significantly larger temperature range 
suggests that we actually observe VRH as predicted by 
AAG Q . Their expression for the magneto- resistivity in 
a 2D superconductor is: 

Mi?(B, T) = ( A)2^„[„^(s)]exp[-(To/T)i/3)] (f) 

where 70 is the vortex conductivity which depends on the 
vortex density n^(B), and Tq is given by: 

T,=K5V{^f (2) 
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where K isa dimensionless factor of order unity, SV is the 
average variation of the pinning energy, is the pairs 
density and Upin is the pinning sites density. It should 
be noted that this theory strictly applies only to the 2D 
case, while the actual film has a quasi 3D nature with 
pancake vortices which have very weak Josephson cou- 
pling in the c direction. The vortices in the film, form 
straight lines which minimize the total pinning potentials 
of the layers. The minimal instanton action (wavefunc- 
tion overlap) for the tunneling events involves motion of 
an individual vortex which has a pinning site in its close 
neighborhood. After the tunneling event, the vortex lines 
up to its new position via classical (above barrier) relax- 
ation. This satisfies the ingredients of the variable range 
hopping model, where the mean tunneling distance in a 
thin film is of the order of I2D / \J ^layers , where I^d is 
the mean distance between uncorrelated pinning sites in 
each 2D layer, and Niayers is the number of layers in the 
film. Therefore, for the quasi 3D case of a thin film, Eq. 
(2) has to be modified by replacing Upm by npi„ x Niayers 
which yields: 

ToifUm) = KSV{ ^ (3) 

^pin-^^ layers 

One can see that for any given field the results of 
Fig. 3 at low temperatures are in good agreement with 
Eq. (1). Ahhough the FTY theory fd\ is applicable 
only at high fields, we also tested their prediction of 
p cx exp [-{To/Tf/^)], but this dependence fit our MR 
data only in the narrow regime of 2-5 K. In the inset 

1 /3 

of Fig. 3, the Tq' coefficient of Eq. (1) is plotted 
versus field. We see that this coefficient is not exactly 
a constant as assumed for low fields by the theory, 
but has a small decreasing contribution of about 30% 
with increasing magnetic field. This behavior can be 
explained by the effect of decreasing barrier height for 
vortex tunneling when the vortex-vortex interaction is 
increased with increasing field. 

Next we make a consistency check of our data at 
H=2 T with Eq. (3). From the data of Fig. 2 just 
below Tc we find that the pinning energy of our films is 
Uo_^ V 550K. For SV of about 10% of V one gets 
6V « 55 K. The measured Tq from the inset of Fig. 3 at 
2 T is 3.8^ w 55 K. Thus, for our data Eq. (3) implies 
that TTUs « npinNiayers- For pairs, rig equals half the 
doping p per copper in the CUO2 plane. In our 60 K 
YBCO phase, p w 0.12 [l3| which yields w 0.06 per 
copper. Thus ripin ~ ttO. 06/170 per copper in a single 
Cu02 plane, where Niayers ~ 170 is the number of Cu02 
planes in the 100 nm thick YBCO film. This implies 
that the average distance between pinning s ites in a 
single Cu02 plane is ho « 3.9/^/71-0.06/170 w 117A 
which is very reasonable, (the 3.9 A here is the in-plane 
lattice constant of YBCO). 

In Fig. 5 we plot the measured magnetoresistance 
MR versus magnetic field for a few temperatures. One 
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FIG. 5: (Color online) Magnetoresistance as a function of 
magnetic field at a few temperatures. The straight lines are 
just connecting the data points, while the parabolic curve is 
a quadratic fit to the data at 10 K for fields of 1-6 T. 

can see that at the low temperatures of 2.5 and 5 K the 
MR is almost linear in field. In contrast, at higher tem- 
peratures the dependence on the field is more complex, 
as seen from the data at 7.5 and 10 K. The linear MR 
dependence on field at a constant temperature can be 
attributed to a constant terminal vortex velocity in the 
YBCO film. This would give rise to an induced voltage 
across the meander line which increases linearly with 
field. The behavior of the MR at higher temperatures 
and higher fields can originate in the weak links which 
are distributed along the meander line. Under these 
conditions, the weak links become normal, and this gives 
rise to an MR which reflects the distribution of the 
strength of the weak links which can be quite complex. 

Finally, we point out that the main difference between 
the present study and the work on the vortex glass 
dynamics [6*], is that here we probed the low field (low 
vortex density) regime. The single vortex tunneling 
theory ,T] predicts a magnetoresistance linear in H 
as observed in Fig. 5, and ignores vortex interaction 
effects which are of higher order in the magnetic field. 
Experimentally, this is justified since at the highest 
field used of 6 T, the distance between vortices is of 
about 20 nm, while the distance between pinning sites 
in the whole film is only 117/ \/170 ~ 0.9 nm. Therefore, 
collective tunneling effects might well be relevant only 
at high fields and temperatures. Since our experiment 
probes the low field regime, (much lower than Hc2), we 
expect single vortex variable range hopping to dominate 
our results. 

In conclusion, vortex variable range hopping in two 
dimensions is consistent with our experimental magne- 
toresistance results in a long meander line of underdoped 
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YBCO thin films at low temperatures. 
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